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PREFACE 


During  the  past  two  and  one-half  years.  The  Rand  Corporation  has 
been  conducting  a  study  for  the  Advanced  Research  Projects  Agency  of 
possible  new  applications  for  lasers.  Several  reports  have  already 
resulted  from  this  study,  which  is  still  in  process.  The  results  pre¬ 
sented  in  these  reports  have  depended,  in  part,  on  detailed  calcula¬ 
tions  of  the  diffraction  patterns  produced  by  focused  iaser  beams. 

These  calculations,  described  herein,  show  relationships  between  char¬ 
acteristics  of  the  beam  as  it  leaves  the  output  aperture  and  those  of 
the  beam  when  it  reaches  a  focus  on  a  distant  surface. 

Some  of  the  general  results  given  in  this  report  have  also  appeared 
elsewhere  in  the  literature.  However,  almost  none  of  our  detailed,  spe¬ 
cific  results  have  been  available  until  now.  It  was  the  need  for  these 
data  in  the  larger  study  that  prompted  us  to  undertake  these  calculations. 
They  are  being  published  here  so  that  they  may  bo  mode  widely  available 
for  use  in  the  design  of  laser  systems  and  in  the  analysis  of  their  per¬ 
formance. 
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SUMMARY 


Some  types  of  lasers  produce  (ideally)  a  monochromatic  collimated 
beam  with  a  gaussian  radial-intensity  variation.  The  beam  can  be  ex¬ 
panded  by  telescope-lika  optics,  projected  from  an  output  aperture, 
focused  onto  a  distant  target,  and  kept  aimed  at  the  same  point  on  the 
target  for  some  appropriate  time  interval.  The  objective  of  this  pro¬ 
cedure  is  to  achieve  a  high  energy  density. 

This  report  examines  two  of  the  factors  that  influence  the  resul¬ 
tant  energy  density:  diffraction,  which  prevents  the  beam  from  focus¬ 
ing  to  a  point,  and  jitter  of  the  diffraction  pattern,  which  spreads 
energies  over  a  larger  area.  Diffraction-pattern  characteristics  are 
found  from  numerical  solutions  of  the  Fraunhofer  diffraction  integral 
for  circular  apertures  with  circular  central  obscurations  and  gaussian 
illuminations.  For  describing  intensity  patterns,  we  have  chosen  the 
"spot"  to  be  the  circular  area  whose  edge  lies  at  the  points  where  the 
Intensity  has  fallen  to  one-half  of  the  peak,  central  value.  Formulas, 
tables,  and  graphs  are  presented  that  enable  vhe  reader  to  estimate  the 
peak  intensity,  the  spot  size,  the  power  within  the  spot,  and  the  av¬ 
erage  intensity  over  the  spot.  Required  inputs  are  the  range,  the 
radiated  power,  the  aperture  diameter,  the  wavelength,  the  fraction  of 
the  aperture  area  that  is  obscured,  and  a  measure  of  the  amount  of 
truncation  of  the  gaussian  illumination. 

The  average  intensities  that  result  from  a  jittering  beam  are 
easily  estimated  if  both  the  jitter  statistics  and  the  diffraction  pat¬ 
tern  are  gaussian — the  two  combine  (convolve)  to  yield  another  gaus¬ 
sian  pattern.  This  report  shows  that  the  high-intensity  portion  of  a 
diffraction  pattern  can  be  fitted  closely  with  a  gaussian  pattern  by 
choosing  the  gaussian  parameters  so  as  to  match  the  size  and  power  of 
the  diffraction  spot. 
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SYMBOLS 


A  = 
D  * 
F  = 

£  3 

f(B)  * 

G  * 

h  * 

I  => 

I  * 
o 

J  a 

O 

m  * 


n  = 
P  - 

R  => 
r  = 
8  = 

t  -• 

u  = 

u  = 
w  = 

o 

z  = 

a  = 


a 


D 


2 

area  within  the  aperture  perimeter  =  irD  /4 
aperture  diameter 

2 

obscuration  fraction  (the  fraction  (e  )  of  area  A  that  is  ob¬ 
scured) 

fraction  of  power  P  in  the  focal  plane  contained  within  the  half¬ 
intensity  circle 

fraction  of  power  P  in  the  focal  plane  contained  within  a  circle 
of  radius  0 

gaussian  power-correction  factor  (GP  is  the  total  power  in  the 
gaussian  approximation  to  a  real  diffraction  pattern) 

dimensionless  angular  radius  of  the  half-intensity  circle  (unit 
angle  =*  X/D) 

intensity  (power  per  unit  area)  at  a  point  in  the  focal  plane 
maximum  intensity  (power  per  unit  area)  in  the  focal  plane 
the  zero-order  Bessel  function 


diffraction  limit  (if  6  is  the  spot  diameter  for  a  diffraction- 
limited  system,  m0  is  the  diameter  for  an  "m  times  diffraction- 
limited"  system) 


aperture  radius  D/2 

taper  parameter  "  -xr - t, -  3  - 

r  r  beam  radius  w 

power  radiated  from  the  aperture  to  the  focal  plane  (no  trans¬ 
mission  losses) 


distance  from  the  aperture  to  the  focal  plane 

radial  distance  in  the  plane  of  the  aperture  (see  Fig.  1) 

dimensionless  standard  deviation  of  a  gaussian  approximation  to 
a  diffraction  pattern  [aD  *  s(X/D)R] 

time  interval  during  which  the  laser  beam  stays  focused  on  a 
target 

wave  amplitude  at  points  in  the  focal  plane  (I  =  luj  ) 

wave  amplitude  at  points  in  the  plane  of  the  aperture 
2 

1/e  radius  of  a  gaussian  beam  incident  on  the  aperture 
2 

1/e  radiu3  of  a  gaussian  beam  at  the  focal  plane  (wq  = 

dimensionless  variable  of  integration  denoting  radial  distance 
in  the  aperture  plane  -  2r/D 

radial  angle  subtended  by  a  radial  distance  e  in  the  focal  plane 
when  viewed  from  the  aperture  =  p/R  (see  Fig  i) 

angular  standard  deviation  associated  with  a  circular  gaussian 
diffraction  pattern  =  o^/R 
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lj  *  angular  standard  deviation  associated  with  a  circular  gausslan 
probability  density  function  =  Oj/R 

0  *  dimensionless  radial  angle  =  (see  Fig.  1  for  a) 

e  =  obscuration  ratio  (the  ratio  of  the  diameter  of  the  central  ob¬ 
scuration  to  the  diameter  of  the  aperture)  =  /f 

A  =  wavelength 

p  »  radial  distance  in  the  focal  plane  (see  Fig.  1) 


a  *  standard  deviation  of  the  two  normal  distributions  that  combiao 
to  form  the  circular  gausslan  distribution 

Op  *  linear  standard  deviation  associated  with  a  circular  gausslan 
diffraction  pattern 

Oj  *  linear  standard  deviation  associated  with  a  circular  gausslan 
probability-density  function  (jitter) 

<j>(8)  *  dimensionless  intensity  at  a  radial  angle  0  in  the  focal  plane; 
unit  intensity  =  PA/A2R2 

d> (S)  ■  dimensionless  intensity  averaged  over  a  disc  of  radius  3  in  the 
focal  plane 

<J>o  *  maximum  dimensionless  intensity  in  the  focal  plane 

<l>  -  dimensionless  intensity  averaged  over  the  half-intensity  disc  in 
the  focal  plane 
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I.  BACKGROUND  AND  OVERVIEW 

Some  types  of  lasers  produce  (ideally)  a  monochromatic  collimated 
beam  with  a  gaussian  radial-intensity  variation. ^  This  beam  can  be 
expanded  by  telescepe-like  optics,  projected  from  an  output  aperture, 
focused  onto  a  distant  target,  and  kept  aimed  at  that  point  on  the  tar- 
pet  for  some  appropriate  time  interval.  The  focusing  produces  a  high 
power  density  by  concentrating  most  of  the  projected  power  onto  a  small 
spot  on  the  target.  Keeping  the  location  of  the  spot  fixed  produces  a 
high  energy  density. 

For  high-quality  optics,  one  of  the  factors  that  determines  how 
much  concentration  of  power  is  possible  is  diffraction.  Section  II  of 
this  report  examines  the  characteristics  of  diffraction  patterns  pro¬ 
duced  by  circular  apertures  with  circular  central  obscurations.  Special 
emphasis  is  placed  on  the  effects  of  truncation  and  obscuration.  A  few 
of  the  patterns  are  described  in  considerable  detail;  i.e.,  intensity 
and  Included  power  are  given  as  a  function  of  radial  distance  in  the 
focal  plane.  For  more  concise  pattern  descriptions  we  have  chosen  the 
"spot"  to  be  the  circular  area  whose  edge  lies  at  the  points  where  the 
intensity  has  fallen  to  one-half  of  the  peak,  central  value.  Formulas, 
tables,  and  grapm*  are  presented  that  enable  the  reader  to  estimate  the 
peak  intensity,  the  spot  size,  the  power  within  the  spot,  and  the  av¬ 
erage  Intensity  over  the  spot  that  can  be  achieved  by  a  given  laser/ 
telescope  system. 

If  the  diffraction  pattern  stays  fixed  with  respect  to  the  target 
surface,  the  resultant  energy  densities  are  easily  determined  by  multi¬ 
plying  the  power  densities  by  the  time  interval.  Section  III  addresses 
the  problem  of  estimating  energy  densities  when  the  center  of  the  dif¬ 
fraction  pattern  moves  or  "jitters'1  about  an  aim  point  on  the  target. 

A  method  for  estimating  the  effect  of  jitter  is  presented  that  depends, 
in  part,  on  approximating  a  real  diffraction  pattern  with  a  gaussian 


pattern.  Adequate  approximations  car  be  constructed  by  using  some  of 
the  pattern  characteristics  developed  in  Section  LI. 

Some  of  tbe  general  results  given  in  this  report  have  also  ap¬ 
peared  elsewhere  in  the  literature.  Mai.v  of  the  rcte.ences  lifted  at 
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II.  DIFFRACTION  PATTERNS 

\ 

This  section  presents  characteristics  of  diffraction  patterns  (or 
"images")  produced  by  circular  apertures  with  circular  central  obscura¬ 
tions.  The  radiation  incident  upon  the  aperture  is  assumed  to  be  a 

monochromatic  beam  with  a  converging  spherical  wave  front  and  a  gauss ian 
* 

radial-intensity  variation.  The  wave  front  leaving  the  aperture  will 
be  truncated  and  sometimes  obscured.  We  will  compute  the  diffraction 
pattern  that  occurs  in  the  geometric  focal  plane,  i.e.,  In  a  plane  that 
is  normal  to  the  beam  axis  and  located  at  the  center  of  crrvature  of 
the  emerging  wave  xront.  Results  of  such  computations  are  bften  called 
"far-field"  results  because  they  also  describe  the  radiation  pattern 
at  great  uxstancaa  from  apertures  illuminated  by  plane  waves.  All  are 

examples  of  Fraunhofer  diffraction.  { 

\ 

INPUTS  AND  METHODOLOGY 

The  emerging  bedm  is  centered  in  the  aperture.  Figure  1  shows 
the  geometry  and  identified  some  of  the  parameters:  D  Is  the  aperture 

\  \ 


Fig.  1 — The  geometry  and  the  coordinqtes 

'i 


throughout  th^s  rep-rt,  intensity  refers  to  power  per  unit  area. 
(In  the  current  literature,  this  quantity  is  often  called  irmdianae .) 
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diameter,  eD  is  the  diameter  o£  the  obscuration,  and  R  is  the  distance 
to  the  focal  plane.  Because  of  the  axial  symmetry  of  the  output  radia¬ 
tion,  the  diffraction  pattern  will  have  the  same  kind  of  symmetry. 

Points  in  the  focal  plane  are  denoted  either  by  their  radial  distance 
p  from  the  center  of  the  pattern  or  by  the  "radial"  angle  a  that  p  sub¬ 
tends  when  viewed  from  the  aperture. 

?oints  in  the  aperture  plane  are  denoted  by  their  radial  distance 
r  from  the  center.  The  gaussian  amplitude  distribution  of  the  aperture 
illumination  can  be  expressed  as 

u(r)  *  uq  exp(-r2/w2)  (1) 

i  \ 

1 

where  w  is  a  constant  defined  to  be  the  beam  "radius."  The  intensity 
2 

is  given  by  |uj  .  Thus  w  is  the  radial  distance  at  which  the  intensity 
2 

has  fallen  to  1/e  of  its  axial  value.  i 

Note  that  Eq.  (1)  describes  the  amplitude  of  a  plane  wave.  If  it 

became  necessary  to  specify  the  wave-front  cur  ature,  u(r)  would  be 

2 

multiplied  by  e  .mplex  phase  factor  of  the  form  exp(-ii»r  /XR).  How¬ 
ever,  Eq.  (1)  is  the  proper  amplitode  expression  to  use  In  the  equations 
that  will  be  presented  later. 


Soluf on  for  an  Infinite  Gaussian  Beam  \ 

An  analytical  solution  exists  for  the  wave  propagation  of  an  ideal 
(infinite)  beam.  ^  Dickson^?  and  Bloom^  have  estimated  that  if  the 
aperture  radius  is,  at  least  twice  as  large  as  the  beam  radius;,  the  ef¬ 
fects  of  the  truncation  are  "negligible"  and  the  analytical  solution 
is  adequate.  For  this  idealized  case,  the  beam  remains  gaussian  every- 

I 

where.  The  diffraction  pattern  at  the  focal  plane  has  a  radius  wo  given 


w 


o 


XR 

TTW 


*The  quantity  t ,  which  is  the  ratio  of  the  two  diameters  eD  and  D, 
is  usually  called  the  obscuration  ratio. 


where  X  la  the  wavelength.  The  focal-plane  intensity  distribution  I(p) 
la  then 


l(p)  -  I  e*p(-2p2/w2) 
O  o 


(2) 


The  intensity  Iq  at  the  colter  of  the  pattern  can  be  found  bv  equating 

the  total  incident  power  P  to  the  integral  of  l(p)  over  the  focal  plane. 

The  resultant  value  of  1  is 

o 


I 

o 


p*w 

2  2 
X  R 


(3) 


The  Diffraction  Equation 

When  the  output  beam  is  significantly  truncated  or  obscured,  the 
diffraction  pattern  is  no  longer  gaussian;  rather,  it  resembles  an  Airy 
pattern,  i.e.,  a  bright  central  disc  surrounded  by  alternate  dark  and 
light  rings.  To  find  the  intensity  at  a  point  in  the  focal  plane,  we 
must  evaluate  the  Fraunhofer  diffraction  integral  which  yields  the  wave 

h  y 

amplitude  U(p).  The  intensity  1  is  then  given  by  1  *  ju!  . 

For  a  uniformly  illuminated  circular  aperture  with  a  central  ob- 
scruation,  the  Fraunhofer  diffraction  integral  is^ 

D 

2 

u(a)  =  Ir  u  f  J0  (nr ttr)  r  dr 
e  2 

where  u  is  the  uniform  wave  amplitude  at  the  aperture,  and  JQ  is  the 
zero-order  Bessel  function.  Note  that  the  integration  is  performed 


To  amplify  an  earlier  comment,  Fraunhofer  diffraction  applies 
either  in  a  focal  plane  or  in  the  far  field  (R  >>  2^/X)  of  a  collimated 
beam. Elsewhere,  the  Fresnel  diffraction  integral  must  be  used. 
Results  for  such  "near-field"  intensities  produced  by  truncated  gaus¬ 
sian  apertures  are  presented  in  Refs,  6  and  7. 
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orer  the  illuminated  portion  of  the  aperture  and  that  the  radial  angle 

a  is  being  used  to  denote  the  radial  location  of  the  point  in  the  dif- 
* 

fraction  pattern. 

If  the  aperture  illuaiuatlon  is  non  uniform,  the  amplitude  u(r) 
must  be  brought  under  the  integral  sign.  For  a  gauss Ian  illumination, 
the  result  is 


D 

2 

0(a)  *  uo  y*  eaqp(-r2/w2)  JQ  (~  or)  r  dr  (4) 

D 

£  2 

The  diffraction  equation  (Eq.  (4))  can  be  cast  in  a  form  that  is 
more  suitable  for  calculating  and  presenting  results.  First,  we  intro¬ 
duce  two  key  aperture  parameters,  n  and  F,  that  specify  the  amount  of 
truncation  and  obscuration.  They  are  defined  as 

*  aPerture  radius  _  D/2 
beam  radius  w 

2 

F  *  the  fraction  of  the  total  circular  aperture  area  (iD  /4)  that 
ir  obscured  *  e2 

(These  parameters  will  be  discussed  in  more  detail  later.)  Next,  we  in- 
troduce  the  dimensionless  radial  znale  S  =  and  change  toe  variable 
of  integration  to  z  *  2r/D.  As  a  result  of  these  four  changes,  the  dif¬ 
fraction  equation  will  become 

>22 

U(6)  »~u  /  enZJ  (rgz)  z  dz 

AK  G  mf  O 

'f 

where  A  is  the  full  circular  area 


The  approxiiratior.  c  *  iR  is  used  here  and  elsewhere  throughout 
this  report. 


Finally,  we  find  an  expression  for  uq  by  equating  the  total  radi¬ 
ated  power  P  to  the  integral  of  the  intensity  over  the  aperture.  That 
is. 


D 

2 


P  -  2  it 


f  % 


o  2  2 

-2n  z  , 
e  r  dr 


*>  9 

.2  /*  -2nV  . 

2Au  I  e  z  dz 

°  J 


/F 


Au  2  0  2 

o  /  -2Fn  -2n  . 
— r  (e  -  e  ) 

2n 


Thus 


so  the  diffraction  equation  becomes 


1 


An  advantage  of  this  formulation  is  that  all  of  the  parameters 

that  have  physical  dimensions  have  disappeared  from  the  integral  and 

2  2 

appear  together  in  the  single  expression  PA/X  R  .  This  quantity  is  of 
special  importance  in  diffraction  theory  because  it  is  the  expression 
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for  the  maximum  intensity  in  the  Fraunhofer  diffraction  pattern  produced 
by  a  uniformly  illuminated  aperture  of  area  A.  It  is  a  natural  choice 
for  a  normalizing  factor,  i.e.,  a  unit  of  intensity.  Therefore,  the 
final  statement  of  the  diffraction  formula  is  as  follows:  We  define  a 
dimensionless  intensity  <}>(0)  by  the  expression 


♦  (B)  - 


2  2 


J“  e  n  Z  Jq  (tt0z)  z  dz 


/F 


(6) 


and  use  <{>(8)  for  describing  the  diffraction-pattern  characteristics. 

The  "physical"  intensity  can  always  be  found  later  from 

1(8)  -  $(B)  •  ~~r  (7) 

XX 

and  linear  radii  in  the  focal  plane  can  always  be  found  from  a  ■  8(X/D) 
and  o  * **  aR. 

The  Maximum  Intensity 

The  maximum  intensity  occurs  at  the  center  of  the  image,  which  cor¬ 
responds  to  8  *  0.  When  8  *  0,  the  function  Jq  ■  1.  This  simplifica¬ 
tion  cf  the  integrand  allows  the  diffraction  Integral  to  be  evaluated 

in  closed  form.  The  resulting  expression  for  the  dimensionless  maximum 
** 

intensity  is 


♦ 


o 


2_  e"Z(1-F)  -  1 
n2  en2(l-F)  +  j 


(8) 


*It  is  also  the  maximum  possible  intensity.  Anv  nonuniformity  of 
the  aperture  illumination  will  reduce  the  maximum. (®) 

**Reference  9  presents  a  closed-form  solution  for  the  on-axis  in¬ 
tensity  for  all  ranges,  including  the  focal  range. 
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£t  is  not  difficult  to  verify  that 


lia  t  -  1  -  F 
n*0  ° 


This  is  a  necessary  result  because  n  0  iaplies  a  uni  for*  aperture  il¬ 
lumination.  In  that  case  we  *ust  have 


o  lV 


where  A*  *  a(1  -  F)  is  the  true  aperture  area. 


luations  for  Included  Power  and  Average  Intensity 


Other  important  diffraction-pattern  characteristics,  aside  fro* 
t(S)  itself,  are  f(3),  the  fraction  of  the  total  power  in  the  focal 
plane  contained  within  a  circle  of  radius  3;  and  s(3),  the  average  in¬ 
tensity  over  a  circle  (or  "disc”)  of  radius  2.  Starting  with 


f(e)«£  f  »ca>  2tc  dt 

J  n 


ana  using  the  relation  c  =  £0/D)R,  we  obtain 


:a)  -f  /  =<a 


)  £  d£ 


To  find  a  useful  expression  for  :(£),  we  start  with  the  required 


relationship 


777T  PA  „  fill  1  p 

* >  72  7 

'“R  '  c“ 
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ann  obtain 


♦  (2) 


2.2 

v  S 


f(3) 


(10) 


Selected  Input  Values  for  n  and  ? 

The  equation  for  f(8)  shows  that  for  a  given  value  of  6,  the  in¬ 
tensity  #  depends  upon  the  values  of  only  two  parameters,  n  and  F.  As 

defined  earlier,  i  *  (0/2) /v.  The  Introduction  of  n  allows  the  aperture- 

2  2 

illumination  amplitude  to  be  written  as  u  *  u  exp [ -n  (2r/D)  ] -  At  the 

°  2 

edge  of  the  aperture  the  intensity  falls  to  exp(-2n  )  of  its  axial  value. 

Figure  2  shows  hew  the  intensity  varies  across  the  aperture  for  n  *  0, 

0.5,  1.0,  1.5,  and  2.0.  These  five  illuminations  encompass  a  very  wide 

* 

range.  Mote  that  n  -  0  implies  a  uniform  illumination.  The  greater 
-he  value  of  n,  the  more  steeply  the  it tensity  curve  tapers  from  the 
center  to  the  edga,  Therefore  we  shall  call  n  the  "taper”  parameter 
in  the  remainder  j£  this  discussion.  This  parameter  is  a  measure  of 
the  non uniformity  of  the  illumination. 

The  obscuration  fraction  F  is  the  fraction  of  the  total  area  A  that 
is  obscured.  Figure  3  shows  the  aperture  appearances  for  F  =  0,  0, 1, 

0.2,  0.3,  0.4,  and  0.5.  Clearly,  they  encompass  a  wide  range  of  obscur¬ 
ations. 

Later  in  this  section  (pp.  13  -  19),  certain  specific  and  aggre¬ 
gate  characteristics  of  diffraction  patterns  will  be  presented  for  all 
combinations  of  n  and  f. 

Proper  Interpretation  of  ?  and  A 

Before  we  present  results,  some  words  of  caution  are  appropriate. 

2  2 

It  is  important  to  remember  that  the  power  P  in  the  expression  PA/X  R 
is  the  power  radiated  from  the  aperture  (to  be  even  more  precise,  P 

i 

By  having  earlier  factored  P  and  D  out  of  the  diffraction  equation, 
we  have  avoided  the  implication  that  no  power  is  radiated  if  n  *  0. 


0.4 


0.6 


0.8 


1.0 


.2 


D/2 


Fig.  2 — Aperture  intensity  distribution 
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F  =  0 


F  =  G.J 


Fig. 3 — Aperture  appearances  for  var’ous  obscuration  fractions,  F 


should  be  the  newer  that  reaches  the  fo<-al  plane,  but  we  will  assume  no 
transmission  losses) .  The  power  of  the  laser  beam  that  enters  the  beam- 
expansion  optics  is  not  necessarily  the  same.  To  formulate  the  diffrac¬ 
tion  equations,  we  postulated  that  the  wave  is  an  ideal  gaussian  before 
it  is  truncated  and  obscured  by  the  aperture.  For  such  a  case,  it  is 
easy  to  show  that  if  the  power  of  the  ideal  beam  were  P',  then  the  power 
transmitted  by  the  aperture  would  be 

P  =  P'  •  (e~2Fn~  -  e“2n  )  (11) 


the  losses  being  caused  by  truncation  and  obscuration.  But  a  real  beam- 
expander  might  not  work  this  way.  This  study  has  not  examined  methods 
for  estimating  such  losses. 


A  second  point  to  reemphasize  is  that  the  area  A  in  the  expression 
PA/X2*2  is  not  necessarily  the  illuminated  area  of  the  aperture.  It  is 
always  equal  to  the  full  circular  area  wD  74.  When  there  Is  obscuration, 
the  illuminated  area  is  equal  to  A(1  -  F). 

RESULTS 

First  we  will  present,  in  detail,  the  patterns  corresponding  to  a 
selected  set  of  four  pairs  of  values  for  n  and  F.  Less  detailed  results 
will  then  be  presented  for  the  30  pairs  considered. 


Detailed  Diffraction-Pattern  Characteristics 

Figures  4  through  6  show  intensity  9  versus  radial  angle  6  for  a 
selected  set  of  tapers  and  obscurations.  The  tapers  are  n  *  0,  corre¬ 
sponding  to  uniform  illumination,  and  n  =  1,  corresponding  to  a  gauss i an 
truncated  at  the  l/ez  radius.  The  obscurations  are  none  (F  =  0),  and  10 
percent  (F  =  0.1).  The  results  in  Figs.  4  thtough  6  were  obtained  by 

direct  numerical  integration  of  the  diffraction  integral.  No  general 

* 

closed-form  solution  exists.  However,  series  solutions  that  may  con¬ 
verge  rapidly  have  been  devised. 

Figure  4  shows  that  tapering  the  illumination  decreases  the  maxi¬ 
mum  intensity,  widens  the  central  disc,  and  reduces  the  intensity  of 
the  rings.  Figures  5  and  6  show  that  obscuration  decreases  the  maximum 
intensity,  narrows  the  central  disc,  and  increases  the  Intensity  of  the 
first  ring. 

Figure  7  displays  the  fraction  of  the  total  power  in  the  image 


*  (5) 

For  uniform  illumination  (n  *  0),  4(B)  can  be  found  analytically. 


The  solution,  first  found  by  Airy,  is 

♦<«>  -  TT-f  l 


2J.  (ir6)  2J  (eirB) 
L  -  F  1 


irB 


cirB 


where  is  the  first-order  Bessel  function.  Results  obtained  from 
this  formula  and  from  the  analytical  formula  for  were  compared  with 
results  obtained  by  numerical  integration  to  assess  the  accuracy  of 
the  numerical  methods.  For  all  comparisons,  quantities  agreed  to  with¬ 
in  at  least  three  significant  figures. 


Intensify,  4> 


Fig  .7 — Fraction  of  image  power  within  disc 

that  lies  within  circles  (i.e. ,  discs)  of  various  radii.  The  curves 
show  the  behavior  of  the  central  portion  of  the  patterr';  they  have  baen 
stopped  at  the  first  null.  If  they  were  extended  farther,  they  would 
resemble  step  functions  as  successive  rings  are  Included  within  the 
circle.  These  curves  illustrate  the  main  effect  of  a  central  obscura¬ 
tion:  It  diffracts  power  into  the  rings. 

Figure  8  displays  the  average  intensity  as  «  function  of  the  size 
of  the  disc  oyer  which  the  average  is  taken.  In  laser-application  stud¬ 
ies,  an  average  intensity  is  often  used  to  describe  the  diffraction 
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Fig.  8 —Average  intensities  v^ithin  disc 

\ 

pattern,  yet  there  Is  no  agreement  among  investigators  about  what  disc 

size  should  be  chosen  for  taking  the  average.  Figure  8  shows  that  the 

average  intensity  can  be  a  strong  function  of  this  choice. 

Table  1  summarizes  some  of  the  major  characteristics  of  the  central 

2 

portions  of  these  four  patterns.  The  full  central  disc,  the  "l/e*-" 

disc,  and  the  "half-intensity"  disc  are  considered.  The  radii  of  the 

latter  two  extend  to  the  points  where  the  intensity  has  fallen,  respec- 
2 

tively,  to  1/e  and  one-half  of  the  peak  value. 
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Table  1 

DI FFIACTIOH-PATTEJBI  CHARACTERISTICS 


Value 

Uniform  illumination 

Gaussian  Illumination 
with  1/e2  Truncation 

Characteristic 

Ho  Obscu¬ 
ration 

102  Obscu¬ 
ration 

Ho  Obscu¬ 
ration 

102  Obscu¬ 
ration 

Peak  intensity9 

0.924 

0.844 

Full  Central  Disc 


Angular  radius*5 

1.22 

1.13 

1.43 

1.20 

Power  fraction 

0.838 

0.668 

0.955 

0.720 

Average  Intensity 

0.229 

0.220 

0.191 

0.203 

1/e2  Disc 


Angular  radius 

0.823 

0.773 

0.915 

0.830 

Power  fraction 

0. 767 

0.615 

0.866 

0.663 

Average  intensity 

0.459 

0.417 

0.419 

0.390 

Half-Intensity  Disc 


Angular  radius 

0.514 

0.488 

0.566 

0.522 

Power  fraction 

0.475 

0.385 

0.530 

0.412 

Average  intensity 

0.727 

0.655 

0.671 

0.613 

^nit  Intensity  *  PA/X2R2. 
^Unit  angle  *  X/D. 


Results  for  the  Half-Intensity  Disc 

The  maxiauu  Intensity  ♦o  is  an  lnportant  descriptor  of  a  diffrac¬ 
tion  pattern,  but  because  it  describes  the  power  density  at  a  single 
point,  it  gives  no  inforaation  about  hew  much  power  can  be  concentrated 
into  some  finite-sized  spot.  For  the  following  presentation  of  results 
we  have  chosen  the  half-intensity  disc  as  the  spot  to  be  described. 

Table  2  shows  the  values  of  the  maximum  intensity  for  all  com¬ 
binations  of  the  taper  parameter  n  and  the  obscuration  function  F.  In 
similar  fashion,  Tables  3  through  5  display  the  radii  of  the  half- intensity 

discs  (denoted  by  h) ,  the  power  fraction  within  them  (denoted  by  f),  and 

—  * 

thier  average  intensities  (denoted  by  $). 

Although  three  quantities  are  presented  for  the  half- intensity 
disc,  only  two  are  independent.  The  three  must  satisfy  the  relation¬ 
ship  < p  -  4f/TT^h2. 
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Table  2 


MAXIMUM  INTENSITY® 


Obscuration 
Fraction,  F 

Maxima*  Intensity,  ♦„ 

Taper 

Parasaater.  n 

0.5 

1.0 

1.5 

2.0 

0.00 

1.000 

0.995 

0.924 

0.719 

0.482 

0.10 

0.900 

0.896 

0.844 

0.682 

0.473 

0.20 

0.800 

0.797 

0.760 

0.637 

0.461 

0.30 

0.700 

0.698 

0.673 

0.584 

0.443 

0.40 

0.600 

0.599 

0.583 

0.523 

0.417 

0.50 

0.500 

0.499 

0.490 

0.453 

0.381 

“unit  Intensity  •  PA/iV. 


Table  3 

HALF- INTENSITY  RADIAL  ANGLE3 


Obscuration 
Fraction,  F 

Half- 

Intensity  Kedlal  Angle,  b 

Taper  Paroaetar.  a 

KSLl 

mam 

KOI 

ESI 

2.0 

0.00 

0.514 

0.526 

0.566 

0.646 

0.773 

0.10 

0.488 

0.496 

0.522 

0.571 

0.639 

0.20 

0.466 

0.471 

0.48J 

0.519 

0.561 

0.30 

0.446 

0.450 

0.461 

0.481 

0.507 

0.40 

0.430 

0.432 

0.439 

0.451 

0.468 

0.50 

0.415 

0.416 

0.421 

0.428 

0.438 

*Unlt  angle  -  1/0. 


Table  4 


FRACTION  OF  IMAGE  POWER  WITHIN  THE 
HALF-INTENSITY  CIRCLE 


Obscuration 
Fraction,  F 

Fraction  of  Iaags  Four 

.  f 

Taper 

Faraawter.  n 

pmag 

tLBm 

msm 

2.0 

0.00 

0.475 

0.494 

0.530 

0.536 

0.514 

0.10 

0.385 

0.396 

0.412 

0.398 

0.347 

0.20 

0.312 

0.318 

0.325 

0.307 

0.260 

0.30 

0.251 

0.254 

0.257 

0.243 

0.205 

0.40 

0.199 

0.201 

0.202 

0.191 

0.164 

0.50 

0.155 

0.155 

0.156 

0.149 

0.132 

Table  5 

AVERAGE  INTENSITY  WITHIN  THE 
HALF-INTENSITY  CIRCLE3 


Obscuration 
Fraction,  F 

Average 

Intaaa 

itr.  J 

Taper. 

Persist 

er.  n 

KflUi 

PQQi 

msm 

msm 

mam 

0.00 

0.727 

0.722 

0.671 

0.521 

0.349 

0.10 

0.655 

0.652 

0.613 

0.495 

0.344 

0.20 

0.582 

0.580 

0.553 

0.463 

0. 335 

0.30 

0.510 

0.509 

0.490 

0.425 

0.322 

0.40 

0.437 

0.436 

0.425 

0.381 

0.304 

0.50 

0.364 

0.364 

0.357 

0.330 

0.278 

*Unlt  lntanalty  •  PA/lV. 
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To  display  the  trends  and  to  facilitate  interpolation,  we  have 

presented  these  results  graphically  in  Fige.  9  through  12. 

As  an  illustration  of  how  these  general  results  can  be  applied  to 

a  specific  case,  consider  a  hypothetical  system  that  radiates  30  kW  of 

continuous-wave  power  at  a  wavelength  of  10.6m  from  an  aperture  30  cm 

in  diameter.  (These  parameters  are  consistent  with  the  performance  of 

CO2  gas-dynamic  lasers  reported  In  Ref.  11.)  We  wish  to  compute  the 

diffraction-pattern  characteristics,  using  the  half-intensity  disc,  when 

the  beam  is  focused  at  a  range  of  500  m.  First,  we  compute  the  unit 
2  2  2 

intensity,  PA/X  R  ,  as  7,550  w/cm  .  In  this  example,  we  assume  that 

2 

the  aperture  illumination  is  gaussian  with  1/e  truncation  and  10  per¬ 
cent  of  the  aperture  area  is  obscured  (n  -  1,  F  «  0.1).  Then,  using 
the  appropriate  numbers  from  either  the  fourth  column  of  Table  1  or 
from  Tables  2  through  4,  we  find 

Peak  intensity  -  0.844  (7,550)  ■  6,370  w/cm^ 

2 

Average  intensity  -  0.613  (7,550)  -  4,630  w/cm 

The  radius  of  the  spot  of  interest  is  0.522  (X/D)R,  which  yields  a  spot 

diameter  of  1.84  cm,  about  the  size  of  a  dime.  The  power  within  the 

spot  is  0.412  P  *  12.36  kW.  It  is  easy  to  verify  that  this  power  di- 

2 

vided  by  the  spot  area  (2.67  cm  )  yields  the  average  intensity  computed 
above. 

It  must  be  remembered  that  the  power  densities  computed  by  this 
method  are  theoretical  upper  limits.  However,  to  illustrate  the  effect 
of  such  a  high  power  density,  if  the  target  were  Indeed  a  silver  dime 
and  if  it  absorbed  10  percent  of  the  incident  radiation,  it  would  melt 
in  less  than  a  second.  (The  reader  is  cautioned  that  melting  silver 
dimes  with  laser  beams  is  neither  cost-effective  nor  legal.) 

An  Empirical  Relationship  Between  Maximum  and  Average  Intensities 

When  the  taper  and  the  obscuration  vary,  the  maximum  intensity 
and  the  average  intensity  f  each  vary  over  a  wide  range.  However,  the 
ratio  ?/$Q  is  very  nearly  a  constant.  We  computed  this  ratio  (7  from 


r 
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Topcr  poranwtar,  n 


Fig.  10 — Half- Intensity  radial  angles 


/X2  R2 )  Froction  of 
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« 


Taper  parameter,  n 


Fig.  11 — Fraction  of  image  power  within  half- intensity  circle 
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Table  5;  from  Table  2)  for  each  of  the  pairs  of  values  for  n  and  F 
and  found  that 


0.  724  s  0.  729  (12) 

*o 

Therefore,  the  relationship  <J>  »  0.727  which  is  correct  for  the  Airy 
pattern  (n  ■  0,  F  -  0) ,  remains  correct  within  0.5  percent  over  the  en¬ 
tire  range  of  tapers  and  obscuration  considered  in  this  study — a  very 
broad  range. 

Apparently,  the  central  portions  of  these  diffraction  patterns — 
out  to  the  half-intensity  circle  at  least — have  a  characteristic  shape. 
That  is,  one  curve  of  $  versus  B  can  probably  be  transformed  very  nearly 
into  any  of  the  others  simply  by  changing  scales. 

The  plausibility  of  this  result  is  indicated  by  the  following  ar¬ 
gument.  The  Bessel  function  J  (x)  can  be  expressed  as  a  power  series 

o  2  ^  ^ 

that  has  the  form  J  (x)  -  1  -  a0x  +  a.x  -  a,x  +  ...»  where  the  a's 

o  i  4  b 

are  known  constants.  If  Jo(ir0z)  in  the  diffraction  integral  (Eq.  (5)) 

were  replaced  by  its  series  and  the  term-by-term  integration  completed, 

the  resultant  expression  for  the  amplitude  U(0)  would  have  the  form 
2  4  6 

U(8)  ■  bg  -  b20  +  b^8  -  b^B  +  . where  the  b's  are  constants.  In 

the  region  of  interest,  8  remains  small;  therefore,  the  first  two  terms 

may  dominate  the  series.  Assuming  this  to  be  true,  it  is  easy  to  show 

2  — 

that  if  U(8)  *  bp  -  b2&  ,  then  ■  0.736,  independent  of  the  values 

of  bQ  and  b2> 

2 

Estimates  for  the  1/e  Disc 

In  light  of  the  pattern  consistencies  described  above,  it  seems 

plausible  that  reasonably  constant  relationships  exist  between  charac- 

2 

teristics  of  the  1/e  disc  and  the  half-intensity  disc.  To  test  this 
conjecture,  we  can  manipulate  the  results  presented  in  Table  1,  which 
shows  characteristics  of  both  discs  for  four  different  sets  of  aperture 


Developed  by  W.  Sollfrey,  of  The  Rand  Corporation. 


-25- 


conditions  .  For  each  set  of  aperture  conditions,  the  ratios  of  results 

2 

for  the  1/e'  disc  to  results  for  the  half-intensity  disc  fall  within 
the  following  narrow  boundaries: 

=*  ratio  of  angular  radii  =  1.60  ±  0.02 
=  ratio  of  power  fractions  *  1.62  ±  0.02 
■  ratio  of  average  intensities  -  0.630  ±  0.007 

Consequently,  the  ratios  do  appear  to  be  reasonably  constant  (but 

the  agreement  between  K.  and  K2  is  probably  a  coincidence).  Therefore, 
2  x  z 

if  the  1/e  disc  is  used  as  the  spot,  its  characteristics  can  be  esti¬ 
mated  from  the  following: 


Angular  radius  ■ 

1.6  h 

(13) 

Power  fraction  ■ 

1.6  f 

(14) 

Average  Intensity  ■ 

0.63  J 

(15) 

where  h,  f,  and  <J>  are  the  values  for  the  half-intensity  disc. 
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III.  ENERGY  DENSITIES  PRODUCED  BY  A  JITTERING  BEAM 

For  most  laser  applications,  it  is  necessary  to  know  how  much  en¬ 
ergy  will  be  delivered  to  certain  areas  of  the  target  when  the  laser 
operates  for  a  time  interval  t.  If  the  diffraction  pattern  stays  fixed 
with  respect  to  the  target  plane,  the  calculation  is  straightforward. 
The  appropriate  intensities  can  be  estimated  by  the  methods  described 
earlier  and  then  multiplied  by  t  to  get  energies.  However,  if  the  pat¬ 
tern  moves  laterally,  the  intensities  at  fixed  points  on  the  target 
plane  do  not  stay  constant  during  t.  This  could  complicate  the  energy 
calculations  by  requiring  Integrations  rather  than  simple  multiplica¬ 
tions. 

In  some  situations,  random  processes  may  cause  the  center  of  the 
pattern  to  "jitter"  about  an  aim  point  during  the  time  interval  t.  The 
problem  of  estimating  the  energy  delivered  by  a  jittering  beam  is  sim¬ 
plified  if  the  time-averaged  Intensity  at  a  point  can  be  treated  as 
being  reasonably  independent  of  the  value  of  t.  It  is  further  simpli¬ 
fied  if  circular  gausslan  functions  can  be  used  to  describe  both  the 
beam  motion  and  the  beam  shape. 

A  METHOD  FOR  ESTIMATING  THE  TIME-AVERAGED  INTENSITY 

» 

Let  the  aim  point  be  the  origin  of  a  rectangular  coordinate  system 
in  the  target  plane.  Let  (x,y)  denote  the  locations  of  the  beam  center 
at  randomly  chosen  instants.  We  postulate  that  x  and  y  are  Independent 
and  are  normally  distributed  with  zero  means  and  Identical  standard 
deviations  Oj.  The  probability  that  the  beam  center  will  be  located 
within  a  small  area  dx  dy  surrounding  a  particular  point  (x,y)  is  thus 
given  by  <j>(x,y)  dx  dy,  where  i|>(x,y)  is  the  product  of  two  normal  proba¬ 
bility  density  functions  (pdf);  i.e. , 

- - —  e-*2/20'  •  j_  e-y2/2°] 

V2"  °j  V2* 


'l'(x.y) 
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or  more  simply. 


<*( P)  -  ~~ J  /2°J  (16) 

2lf0j 

2  2  2 

where  p  ■  x  +  y  .  The  function  ♦(p)  is  called  the  circular  gaussian 
pdf.  For  want  of  a  better  term,  Oj  is  often  called  the  "standard  devi¬ 
ation"  of  the  circular  distribution,  even  though  it  is  not  the  standard 
deviation  of  p. 

It  was  shown  in  Eqs.  (2)  and  (3)  that  for  a  gaussian  beam,  the  in¬ 
tensity  I(q)  at  a  radial  distance  q  from  the  axis  of  the  beam  can  be 
expressed  in  the  form 


I(q) 


o  2.  2 

_  -2q  /w 

I  e  o 

o 


2 

where  I  ■  2  P/irw  ,  and  v  is  the  beam  "radius"  at  the  focus.  For  con- 
o  o  o 

sistency  of  notation,  it  is  advantageous  to  introduce  a  beam  '^standard 
deviation"  by  the  substitution  oQ  *  w^/2.  The  resulting  expression 
for  I(q)  is 


I(q) 


(17) 


which  is  identical  in  form  to  a  circular  gaussian  pdf  except  for  the 
inclusion  of  the  power  P. 

As  the  beam  moves,  the  intensity  at  a  fixed  point  (xQ,yo)  will 
vary.  When  the  center  of  the  beam  is  at  a  point  (x,y)  the  Intensity 

« <vV  1* 


I(x0,y0.x,y)  - 


2w  at 


•*P  <  - 


(X  -  x)J  + 
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To  find  the  expected-value  or  "average”  intensity  I(xQ,yo)  at  the  point 

{x  ,y  }  we  must  evaluate  the  integral 
o  o 


00  oc 

l(xo,yQ)  »  //  I(x0,j0,x,y)  *  <Kx»y)  dx  dy 

oc  oc 


In  Ref.  12  it  is  shown  that  this  integral  yields  the  result 


+  °j'> 


e-%'2(4  + 


0J> 


(18) 


where  pQ  *  +  y^.  This  is  another  gaussian  distribution.  It  indicates 

that  the  formula  for  the  average  intensities  has  the  same  mathematical 
form  as  it  would  have  if  the  gaussian  diffraction  pattern  had  been  held 
stationary  but  broadened  so  that  its  "standard  deviation"  becamevo''  +  o^. 
This  is  a  most  useful  result  for  two  reasons:  (1)  The  energy  calcula¬ 
tions  reduce  again  to  simple  multiplications  by  t,  and  (2)  gaussian 
functions  are  usually  easy  to  manipulate.  For  example,  for  a  gaussian 

intensity  distribution  with  total  power  P  and  standard  deviation  a,  it 

2 

is  straightforward  to  show  that  the  maximum  intensity  is  I  ■  P/2tto  , 

the  radius  of  the  1/e^  disc  is  2o,  the  power  within  it  is  0.865  P,  and 

02) 

its  average  intensity  is  0.432  I  .  For  the  half-intensity  disc, 

o 

the  corresponding  quantities  are:  radius  *  1.177  a,  contained  power 
«  0.5  P,  and  average  intensity  =  0.721  Iq. 


THE  GAUSSIAN  APPROXIMATION 

To  use  the  simplified  treatment  of  jitter  described  above,  we  must 

•k 

approximate  a  real  diffraction  nattern  with  a  gaussian  pattern.  Clearly, 
the  central  portions  where  the  high  intensities  occur  are  the  regions  we 
should  try  to  match  most  closely.  It  follows  that  a  "good"  gaussian 


"Real"  means  only  that  a  pattern  has  been  computed  from  a  diffrac¬ 
tion  integral. 
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approximation  will  not  necessarily  represent  the  delivery  of  the  same 
total  power  as  is  contained  in  the  real  pattern.  Therefore,  we  will 
reserve  the  symbol  P  for  the  real  total  power  and  define  a  "gaussian 
power-correction  factor"  G  such  that  GP  is  the  total  power  in  the  gaus¬ 
sian  image.  The  formula  for  the  gaussian  intensity  pattern  t'  en  becomes 


I(P) 


GP 

2it0D 


e 


-P 


2 


(19) 


A  more  convenient  form  for  this  equation  can  be  achieved  by  intro¬ 
ducing  the  same  kinds  of  dimensionless  quantities  as  were  used  in  Sec¬ 
tion  II.  The  variable  8  and  a  "dimensionless  standard  deviation"  s 
are  introduced  by  using  the  relationships  p  =  8()/D)R  and  rr^  «=  s(\/D)R. 

The  resulting  formula  for  the  intensity  is 


1(8)  = 


PA 

2  2 
X  R 


2  2 
2G  -6 72s 

2  2  C 
IT  S 


(201 


2  2 

Then,  by  again  choosing  the  unit  of  intensity  to  be  PA/X  R  ,  we  obtain 
the  dimensionless  gaussian  intensity  distribution  $(d) 


Me)  = 


9  2 

2G  -B"/2s 

2  2  6 
TT  8 


(21) 


In  this  formulation,  all  of  the  characteristics  ot  the  gaussian 
pattern  are  determined  by  the  values  of  the  two  independent  parameters 
s  and  G.  By  choosing  an  appropriate  pair  of  values  for  these  parameters, 
we  car  match  at  least  two  of  the  characieiisti cs  of  a  real  pattern. 


Methods  for  Selecting  s  and  G 


Section  II  presented  characterise  ii-«  ot  real  dirtra  *ilon  patterns 


-30- 


that  result  from  various  tapers  and  obscurations.  In  that  discussion, 

we  emphasized  the  maximum  intensity  and  three  characteristics  of  the 

half-intensity  disc:  the  radial  angle  h,  the  power  fraction  f,  and  the 

average  intensity  4> .  Also  it  was  noted  that  the  three  half-intensity 

-  22 

quantities  are  not  independent  because  of  the  relationship  1  *  If  h 
Consequently,  4>0»  h,  f,  and  <i>  represent  only  three  irdependent  charac¬ 
teristics. 

A  gaussian  approximation  can  fit  two  of  these  three  independent 
characteristics.  The  method  consists  of  first  using  Eq.  (21)  to  express 
<|>o,  h,  f,  and  $  for  a  gaussian  in  terms  of  s  and  G.  The  maximum  inten¬ 
sity  is 


*  .  2€_ 

^o  2  2 
it  s 


(22) 


To  find  the  half-intensity  radial  angle  h,  we  set  exp(-bZ,2 s^j  = 
and  obtain 


h  *  8  •  /in  4  *  1.177  s  (23) 

The  fraction  f  of  the  real  image  power  P  witnin  the  half-intent ifv 
circle  Is 


/■ 


<K6)  p  d6 


•  2/o  2 
G(1  -  e_a  /2S  l 


-  j 


Finally,  *  he  average  intensity  f  is 


i  ■- 


2  > 
it 
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jr 


♦  «  - - -  0.721  (26) 

ir  s  »,n  4  ns 

Tbe  above  relationships  can  be  used  in  reverse  to  select  s  and  G 
when  h,  f,  and  $  are  given  for  a  real  pattern.  To  make  the  gauss,* an 
pattern  have  the  same  half- intensity  radius  h  as  the  real  pattern,  we 
set 


h 

S  *  1.177 


(27) 


The  parameter  G  can  then  be  chosen  so  that  either  <J>0  or  <f>  is  matched. 
To  match  the  peak  intensity, 


G  = 


(28) 


If  we  wish  to  match  the  average  intensity  <j>,  the  proper  formula  for  G 
is  either 


G  = 


2  2 
IT  S 


in  4  - 

-  $ 


(29) 


or 


G  =  2f 


(30) 


whichevei  i3  more  convenient. 

A  Numerical  Example  and  Recommended  Method  , 

Consider  the  diffraction  pattern  that  results  from  the  aperture 

conditions  r.  *  l  and  F  *  0.1.  (This  is  a  gaussian  illumination  with 
2 

1/e  truncation  and  10  patient  obscuration,  which  was  examined  in  de¬ 
tail  in  Section  It.)  ^rom  'table  1  (or  Tables  2  through  5),  we  find 
$  =«  0.844,  h  -  0,522,  f  =  0.412,  and  q>  »  0.613. 
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We  approximate  It  with  a  gaussian  that  has  the  correct  radius  and 
average  intensity  fot  the  half-intensity  disc;  i.e.. 


h 

sa>  17177 


0.444 


C  -  2f  -  0.824 


For  these  values  the  equation  for  the  gaussian  intensity  distribution 
(Eq.  (21))  becomes 


♦  (6) 


0.847 


-B2/0.394 

e 


(31) 


and  the  resultant  fraction  f(g)  of  the  real  image  power  P  that  is  with¬ 
in  a  circle  of  radius  8  can  easily  be  shown  to  be 

f (8)  -  G ( 1  -  /2s  ) 


=  0.824 


(1 


-3  /0. 394. 
e  ) 


(32) 


In  Figs.  13  ar.d  14  the  intensities  and  the  power  fractions  of  the 
approximation  are  compared  with  those  of  the  real  diffraction  pattern. 
This  comparison  shows  that  the  gau3siar.  function  is  an  excellent  ap¬ 
proximation  for  the  central,  hi^h-intensity  portions  of  the  pattern. 

Note  that  although  we  made  no  explicit  attempt  to  match  the  peak 
value  4'0 ,  a  nearly  perfect  match  occurred  (0.844  for  the  real  pattern; 
0.847  for  the  gaussian  pattern).  This  agreement  is  not  a  coincidence, 
and  it  is  not  limited  to  the  specific  example  chosen  here.  In  Section 
II  (Eq.  (12)),  it  was  found  that  for  real  diffracrior.  patterns,  the 
relationship  <p  -  0.727  $  remains  correct  to  within  about.  O.b  percent. 
For  a  gaussian  pattern,  it  is  easily  seen  from  Fqs,  (22)  and  (26)  that 

<t  =  0.721  1  ,  These  two  equations  are  similar  enough  chat  it  makes 
o  o 

little  difference  whether  the  gaussian  is  fitted  to  ‘  ^  t  . 


Flection  ot  Jnxxre  power  w*thin  half- intensity  circle 
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Fig.:  13 — Comparison  of  intensities 


t'ici,  14 — Comparison  of  power  fractions 
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The  simplest  approach  appears  to  be  the  one  adopted  for  this  ex¬ 
ample.  Namely,  choose  s  and  G  so  that  the  gaussian  pattern  and  the 
real  pattern  yield  the  same  radius  and  average  intensity  for  the  half¬ 
intensity  disc.  The  formulas  are 


s 


h 

1.177 


G  =■  2f 

where  h  and  f  are  as  given  in  Section  II. 


Contents  on  a  Special  Case 

An  aperture  illumination  that  is  often  chosen  as  "representative" 

2 

is  an  unobscured  gaussian  with  1/e  truncation  (n  *  1,  F  **  0).  The 
resultant  diffraction-pattern  characteristics  for  this  case  are  pre¬ 
sented  in  the  third  column  of  Table  1  in  Section  II. 

Suppose  we  wish  to  approximate  this  pattern  with  a  gaussian  that 
2 

matches  the  1/e  disc  rather  than  the  half-intensity  disc.  From  Table  1, 
the  disc  radius  (2s)  is  0.915,  so  we  let  s  ■  0.915/2  -  0.458.  The  power 
fraction,  from  Table  1,  is  0.866.  By  coincidence,  this  is  almost  ex¬ 
actly  the  correct  value  for  a  gaussian  (0.865),  and  therefore,  no  power 
correction  is  needed  (G  -  1).  Consequently,  a  gaussian  distribution 
with  a  total  power  P  and  a  standard  deviation  of  c?D  -  0.458  —  R  will 
match  the  radius  and  the  average  intensity  of  the  1/e^  disc  of  the  real 
pattern. 

Now  consider  a  totally  different  derivation  of  a  gaussian  pattern. 

An  approximation  that  has  sometimes  appeared  in  the  literature  results 
from  matching  P  and  choosing  so  that  the  gaussian  has  the  same  maxi¬ 
mum  intensity  as  that  of  the  ideal  Airy  pattern.  Equating  the  two  maxima 
gives 


ir 
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P  D 

xV 
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which  yields 


o 


D 


(33a) 


-  0.450  (33b) 

which,  by  another  coincidence,  is  very  close  to  *•*„■»  value  of  o^  derived 
from  Table  1. 

These  coincidences  help  to  explain  and  to  Justify  equations  devel¬ 
oped  by  Peckham  in  Ref.  12.  In  that  report  he  postulates  the  same  ap- 

2 

erture  illumination:  gaussian  with  1/e  truncation.  He  then  Introduces 

Eq.  (33a),  attributes  it  to  Buck  (Ref.  13),  and  states  that  it  predicts 
2 

the  1/e  radius  of  the  real  diffraction  pattern.  Finally,  he  conjec- 

2 

tures  that  a  gaussian  with  the  same  1/e  radius  (and  the  same  power) 

* 

will  be  an  adequate  approximation  to  the  real  image. 

The  analyses  presented  here  show  that  a  gaussian  with  the  same 
2 

power  and  1/e  radius  as  the  real  image  would  also  have  the  same  aver- 

2 

age  intensity  over  the  1/e  disc.  Therefore,  they  support  Peckham' s 
conjecture  that  this  approximation  is  adequate  (for  this  particular 
aperture  illumination) . 

The  method  recommended  in  the  present  report  produces  a  slightly 
different  approximation.  The  power  in  the  gaussian  becomes  1.06  P  and 
the  standard  deviation  increases  to  ■  0.481  R.  The  high- intensity 
portion  of  the  pattern  is  fitted  more  closely,  but  it  seems  doubtful 
that  these  differences  would  prove  to  be  Important. 

FINAL  FORM  OF  THE  INTENSITY  EQUATION 

For  the  gaussian  intensity  distribution  that  results  from  a  jitter¬ 
ing  beam,  the  average  physical  Intensity  over  the  spot  is  kIQ,  where  IQ 
is  Che  maximum  intensity  and  k  is  a  constant  that  depends  only  on  the 


In  subsequent  personal  communications,  Peckham  explained  that  he 
was  aware  of  Eq.  (33a)  and  discovered  that  it  was  consistent  with  numer¬ 
ical  results  presented  graphically  by  Buck,  but  ha  was  unaware  of  its 
derivation. 
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definition  of  spot  size.  As  noted  earlier,  k  -  0.432  for  the  1/e2  disc 
and  0.721  for  the  half-intensity  disc.  Here  we  need  only  discuss  IQ. 
From  Eq.  (18),  but  with  P  replaced  by  GP,  the  proper  formula  for 

I  is 
o 

I - GP _ 

°  21r(o2  +  c2) 


In  this  formulation,  aD  and  Oj  have  the  dimensions  of  length.  Since 
jitter  and  the  beamspread  are  usually  treated  as  angles  viewed  from  the 
aperture,  we  introduce  the  angular  standard  deviations  and  ctj  to  get 


I 

o 


GP _ 

2uR2 (a2  +  «$) 


and  since 


the  equation  becomes 


(34) 


I 


o 


_ GP _ 

2*R2[(a£)2  +  ajl 


(35) 


This  is  a  basic  equation  for  the  intensity.  It  explicitly  contains  the 
parameters  G  and  s  that  are  chosen  to  fit  the  nonjittering  pattern  with 
a  gauss lan.  However,  one  alteration  is  needed  to  allow  the  analyst  to 
introduce  a  degradation  in  beam  quality:  The  methods  used  to  obtain  the 
results  presented  in  Section  II  imply  that  the  beam  is  "diffraction- 
limited."  In  Ref.  12  an  "m  times  diffraction-limited"  beam  is  defined 
as  one  for  which  the  standard  deviation  of  the  beam  is  Increased  from 
dp  to  mo^.  Thus,  introduction  of  this  parameter  gives  the  final  form 
of  the  intensity  equation, 


I 

o 


_ GP 

,  _2,  2,  X.2 

2ttR  [m  (s  -) 


(36) 


(12) 

This  equation  is  similar  in  form  to  one  developed  by  Peckham' 
and  used  by  Davis to  study  the  relative  advantages  of  laser  systems 
that  produce  different  wavelengths. 
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Equation  (36)  can  be  written  in  many  different  ways,  some  of  which 
might  not  be  immediately  recognizable  as  the  same  equation.  For  example, 
R.  0.  Hundley  of  Rand,  in  several  publications  prepared  for  ARPA,  ex¬ 
presses  the  intensity  as  an  ideal  value  multiplied  by  a  degradation 
factor  that  accounts  for  jitter  (and  degradation  associated  with  at¬ 
mospheric  effects).  To  convert  to  Hundley’s  equation,  we  let  m  ■  1 

and  rewrite  I  ■  k  I  as 

o 


2  2  2  2 

where  K  =  kG/2ns  .  The  term  KPD  /X  R  is  the  undegraded  intensity.  If 
I  is  to  be  the  average  intensity  within  the  half-intensity  circle,  then 
K  is  simply  J  tt/4,  where  is  found  from  Section  II.  One  other  change 
is  needed  to  complete  the  conversion.  As  long  as  the  correct  value  for 
the  ratio  cXj/a^  is  maintained,  other  measures  of  the  jitter  statistics 
or  beam  dimensions  can  be  used.  Hundley  uses  0j,  the  angular  diameter 
of  the  circle  that  contains  the  beam  center  50  percent  of  the  time,  and 
0D,  the  angular  diameter  of  the  half-intensity  circle.  These  are  con¬ 
sistent  parameters.  He  notes  that  0_  «*•  X/D  for  a  wide  range  of  aperture 

^  * 

conditions  (this  corresponds  to  h  *  0.5)  and  uses  0j/(X/D)  for  the  ap¬ 
propriate  ratio.  The  resulting  equation  is 


I  -  K 


P  D 

2  2 
\  R 


1  + 


(X/D)' 


* 


See  Fig. 


10. 
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